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Abstract 

Graphene holds great promise for post-siUcon 
electronics, however, it faces two main challenges: 
opening up a bandgap and finding a suitable sub- 
strate material. In principle, graphene on hexago- 
nal boron nitride (hBN) substrate provides poten- 
tial system to overcome these challenges. Recent 
theoretical and experimental studies have pro- 
vided conflicting results: while theoretical stud- 
ies suggested a possibility of a finite bandgap of 
graphene on hBN, recent experimental studies 
find no bandgap. Using the first-principles density 
functional method and the many-body perturba- 
tion theory, we have studied graphene on hBN 
substrate. A Bernal stacked graphene on hBN has 
a bandgap on the order of 0.1 eV, which disap- 
pears when graphene is misaligned with respect 
to hBN. The latter is the likely scenario in real- 
istic devices. In contrast, if graphene supported 
on hBN is hydrogenated, the resulting system 
(graphone) exhibits bandgaps larger than 2.5 eV. 
While the bandgap opening in graphene/hBN is 
due to symmetry breaking and is vulnerable to 
slight perturbation such as misalignment, the gra- 
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phone bandgap is due to chemical functionaliza- 
tion and is robust in the presence of misalignment. 
The bandgap of graphone reduces by about 1 eV 
when it is supported on hBN due to the polariza- 
tion effects at the graphone/hBN interface. The 
band offsets at graphone/hBN interface indicate 
that hBN can be used not only as a substrate but 
also as a dielectric in the field effect devices em- 
ploying graphone as a channel material. Our study 
could open up new way of bandgap engineering in 
graphene based nanostructures. 
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Graphene exhibits remarkable electronic proper- 
ties compared to the conventional materials such 
as Si and III-Vs making it an attractive material 
for next generation electronic devices. ^ For practi- 
cle applications, devices made from an atomically 
thin material such as graphene should be supported 
on a substrate. Typically graphene devices are fab- 
ricated on Si02 substrate, however, carrier mobil- 
ity in graphene on Si02 reduces due to charged 
surface states, surface roughness, and surface opti- 
cal phonons in Si02.^'^ Several other oxide-based 
substrates have been investigated so far, how- 
ever, none yields significant improvement over 
Si02.^ Recently, graphene supported on a hexago- 
nal boron nitride (hBN) substrate was found to ex- 
hibit much higher mobility compared to any other 
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substrate. ^'^ High mobility of graphene on hBN is 
enabled by extremely flat surface of hBN and ab- 
scence of dangling bonds at the graphene/hBN in- 
terface. ^ 

Another important challenge for the use of 
graphene in devices is the lack of controllable 
bandgap. ^ A bandgap can be opened through 
quantum confinement by patterning graphene into 
the so-called graphene nano-ribbons (GNRs).^'^ 
However, it is difficult to control the bandgap 
in GNRs due to its sensitivity to the width and 
edge geometry.^ Alternatively, the bandgap can be 
opened by chemical functionalization of graphene 
with a variety of species such as H, F, OH, 
etc.^ The hybridization of the functionalized C 
atom changes from sp^ to sp^, which opens up 
a bandgap. The bandgap opened by this mecha- 
nism is expected to be more robust in the presence 
of disorder compared to the bandgap opened by 
the quantum confinement. ^ The bandgap opening 
by H-functionalization/hydrogenation of graphene 
has been a subject of several recent experimen- 
tal and theoretical studies, which show that the 
bandgap of graphene can be tuned by controlling 
the degree of hydrogenation. ^^"^^ 

Here we report the electronic structure of 
graphene and single-sided hydrogenated graphene 
(graphone) supported on the hBN substrate calcu- 
lated using the first-principles density functional 
method and the many-body perturbation theory 
in the GW approximation, the state-of-the-art 
method for accurate predictions of the electronic 
structure. Theoretical studies have suggested a 
possibility of inducing a bandgap in graphene 
when supported on the hBN substrate, how- 
ever, recent experimental studies find no bandgap 
in this system. ^'^ Earlier calculations based on 
the tight-binding model ascribe this discrepancy 
to the random stacking arrangement of graphene 
on hBN.^ The tight-binding model can not be 
reliably used for quantitative predictions involv- 
ing novel materials such as graphene on hBN 
especially because the interlater hopping param- 
eters between graphene and hBN are not known. 
First-principles calculations are therefore required 
for accurate prediction of electronic structure of 
such novel material systems. Our calculations 
show that slight misalignment of graphene closes 
the bandgap induced by hBN in Bernal stacked 



graphene. We consider hydrogenation as an alter- 
native to induce bandgap in graphene supported 
on hBN. Hydrogenation of a substrate supported 
graphene results in 50% hydrogen coverage and 
the resulting material is called as graphone. ^'^^ 
Free standing graphone exhibits a bandgap larger 
than 2.5 eV. Interestingly, the bandgap of gra- 
phone reduces by about 1 eV due to the substrate 
induced polarization effects ^^"^^ when it is sup- 
ported on hBN. However, unlike graphene, the 
bandgap of graphone is unaffected by the mis- 
alignment with respect to hBN. The band offsets 
at graphone/hBN interface suggest that hBN can 
be used not only as a substrate but also as a dielec- 
tric in the field effect devices employing graphone 
as a channel material. 

The electronic structure calculations are per- 
formed in the framework of density functional the- 
ory (DFT) within the local density approximation 
(LDA) as implemented in the ABINIT code.^^ The 
Trouiller-Martins norm-conserving pseudopoten- 
tials^^ and the Teter-Pade parameterization for the 
exchange-correlation functional are used. To en- 
sure negligible interaction between periodic im- 
ages, a large value (10 A) of the vacuum re- 
gion is used. The Brillouin zone is sampled us- 
ing Monkhorst-Pack meshes of different size de- 
pending on the size of the unit cell: 18 x 18 x 1 
for Bernal stacked graphene on hBN, 6x6x1 
for misaligned graphene on hBN, 18 x 18 x 1 for 
chair-graphone, and 8 x 8 x 1 for boat-graphone. 
For the plane wave expansion of the wavefunc- 
tion, a 30 Ha kinetic energy cut-off is used. The 
quasiparticle corrections to the LDA bandstructure 
are calculated within the GqWo approximation and 
the screening is calculated using the plasmon-pole 
model. 

The Vienna ab initio simulation package 
(VASP),^^ which provides well-tested implemen- 
tation of van der Waals interactions (vdW),^^ 
is used to compute the equilibrium distance be- 
tween graphene (or graphone) and hBN. The PAW 
pseudopotentials,^^ the PBE exchange-correlation 
functional in the GGA approximation,^^ and the 
DFT-D2 method of Grimme^^ are used. Same 
values of energy cutoff, vacuum region, and k- 
point grid as in the ABINIT calculations are used 
in VASP calculations. The optimized geometries 
calculated using VASP and ABINIT without in- 
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Figure 1 : (a) The LDA and GW bandstructures of 
Bernal stacked graphene on hBN substrate. The 
inset shows a small bandgap opening at the K 
point, (b) Bandgap closing due to misalignment 
illusrated by the LDA density of states and GW 
bandgaps of the heterogeneous bilayer of graphene 
and hBN in perfect Bernal and misaligned stack- 
ing arrangements. The inset shows atomistic 
schematic of a commensurate unit cell of the mis- 
aligned bilayer. 

eluding vdW interactions are found to agree well 
with each other. 

Figure [figure] [1][] 1(a) shows LDA and GW 
bandstructures of Bernal stacked graphene on hBN 
substrate such that half of the C atoms in graphene 
are positioned exactly above the B atoms. This 
stacking arrangement has been found to be low- 
est energy configuration.^^ The equlibrium dis- 
tance between graphene and hBN is 3. 14 A. Three 
monolayers of the semi-infinite hBN substrate are 
included in the simulation domain to ensure that 
the GW bandgap is converged. The bands con- 
tributed by hBN are identified separately. The 
weak interlayer interaction between graphene and 
hBN allows graphene to retain its linear bandstruc- 
ture near the K point. Underlying hBN substrate 
induces sublattice asymmetry on the graphene lat- 
tice opening up a small bandgap at the K point as 
shown in the inset of Figure [figure] [1][] 1(a). The 
LDA bandgap is smaller by about 0.1 eV com- 
pared to the GW corrected bandgap. 

Contrary to the above conclusion and earlier the- 
oretical studies, recent transport measurements on 



graphene supported on hBN show no evidence 
of bandgap. ^'^ The perfect Bernal stacking of 
graphene on hBN is difficult to achieve in the ex- 
periments and random orientation is more proba- 
ble. To mimic the random orientation, we simu- 
late larger supercells where the stacking between 
the graphene layer and the underlying hBN sub- 
strate deviates from the ideal Bernal stacking. To 
reduce the computational requirements only one 
layer of hBN is included. In the Bernal stacking, 
the graphene layer is rotated by an angle = 30° 
with respect to hBN. We consider three rotation 
angles 21.8°, 32.2°, and 13.2°. The commensu- 
rate supercells for these rotations contain 28 (14 
C, 7 B, 7 N), 52 (26 C, 13 B, 13 N), and 76 (38 
C, 19 B, 19 N) atoms respectively. The super- 
cell for 21.8° rotation is depicted in the inset of 
Figure [figure] [1][] 1(b). The remaining two su- 
percells and the Brillouin zones of all three su- 
percells are shown in the supplementary material. 
The commensuration conditions derived in^^ are 
used to generate the supercells. In these supercells, 
the sublattice asymmetry induced by the hBN sub- 
strate on the graphene layer is significantly re- 
duced compared to the Bernal stacking. 

The LDA density of states of the misaligned 
graphene {0 = 21.8°) is compared with that 
of the Bernal stacked graphene in Figure [fig- 
ure] [1][] 1(b). The Bernal stacked graphene has 
a finite bandgap, which closes due to the slight 
misalignment. Since LDA is known to under- 
estimate the bandgap, we have calculated the 
GW corrections to the bandgap of misaligned 
graphene. The GW corrected bandgap of the 
Bernal stacked graphene increases from E^^^ — 
68 meV to E^"^ = 145 meV while the GW cor- 
rected bandgap of the misaligned graphene re- 
mains 0. The GW bandgaps for rotations 32.2° 
and 13.2° are also 0. The Moire patterns larger 
than those shown in Figure [figure] [1][] 1(b) were 
recently observed in graphene supported on hBN.^ 
These graphene samples indeed showed bandgap 
similar to our calculations. 

The above calculations indicate that due to mis- 
alignment hBN substrate can not reliably induce 
a bandgap in graphene and alternative approaches 
are required to open up a sizable bandgap. Recent 
experimental measurements and theoretical calcu- 
lations show that the bandgaps on the order of 1 eV 
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Figure 2: Atomistic schematics and bandstructures 
of graphone in (a) chair and (b) boat conforma- 
tions. Unit cells are depicted by the dotted lines in 
the top views of atomistic schematics while Bril- 
louin zones are shown in the insets of bandstruc- 
ture plots. 

can be opened up in hydrogenated graphene. ^^"^^ 
Furthermore, the bandgap can be tuned by con- 
trolling the degree of hydrogenation. Now, the 
question arises, what is the effect of hBN sub- 
strate on the electronic structure of hydrogenated- 
graphene? To address this issue, we consider 
single- sided semi-hydrogenated graphene, which 
is referred to as graphone. Recent experimental 
and theoretical study has shown that hydrogena- 
tion of a substrate supported graphene results in 
50% hydrogen coverage, which we have used 
for the present study. Graphone has two distinct 
configurations chair (Figure [figure] [2] [] 2(a)) and 
boat (Figure [figure] [2] [] 2(b)), depending on the 
relative placement of hydrogenated C atoms. Both 
chair and boat configurations are energetically sta- 
1^1^8,13,16 ^j^^y synthesized by differ- 

ent hydrogenation processes. A direct single-sided 
hydrogenation of graphene is likely to result in 
boat-graphone, which is more stable compared 
to chair-graphone.^'^^ On the other hand, chair- 
graphone can be synthesized by selectively des- 



orbing hydrogen from one side of graphane.^^'-^^ 
The hybridization of C atom changes from sp^ to 
sp^ upon hydrogenation and the planar structure of 
graphene becomes non-planar as depicted in Fig- 
ure [figure] [2] []2. 

Graphene is weakely bonded to hBN by van 
der Waals interaction and the frontier orbitals of 
graphene, which take part in bonding with hydro- 
gen are virtually unaffected by hBN. Therefore 
the stable configurations of graphone are likely 
to be unaffected by hBN. Indeed, the optimized 
atomic structures of graphone on hBN are virtu- 
ally identical to earlier studies, ^'^^'^^ which did 
not include any substrate. The situation is differ- 
ent with the substrates, which have stronger inter- 
action with graphene. For example, graphene on 
Ir(lll) substrate exhibits specific hydrogenation 
patterns based on the local orientation of graphene 
with respect to Ir(l 1 1) surface. Hydrogen atoms 
may desorb or diffuse through graphone into the 
hBN substrate at elevated temperatures thereby af- 
fecting the local electronic struture of graphone. 
However, such phenomena require in depth analy- 
sis using methods such as the molecular dynamics 
and are out of the scope of the present work. 

Chair and boat-graphone exhibit very different 
electronic and magnetic properties as evident from 
their bandstructures (Figure [figure] [2] [] 2). The 
ground state of chair-graphone is ferromagnetic 
while boat-graphone has a nonmagnetic ground 
state. The bandstructure of chair-graphone is 
highly spin polarized with a spin splitting of 
2.79 eV, which is also its bandgap. This fea- 
ture makes chair-graphone an attractive material 
for spintronics.^ Boat-graphone has a spin de- 
generate bandstructure and behaves as an insu- 
lator with a large bandgap of 5.15 eV. Hydro- 
gen atoms may migrate over the graphone lattice 
resulting in the combination of chair- and boat- 
configurations. In such structures magnetic order- 
ing of chair-graphone may not be preserved. The 
bandgap is, however, expected to persist.^ 

The electronic structure of an atomically thin 
material such as graphone, when it is used in de- 
vices, is expected to be highly dependent on the 
surrounding materials. The atomistic schemat- 
ics of graphone supported on hBN substrate and 
calculated bandstructures of graphone-hBN super- 
cells are shown in Figure [figure] [3] [] 3. The gra- 
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Figure 3: Atomistic schematics and bandstructures 
of (a) chair-graphone and (b) boat-graphone sup- 
ported on hBN substrate. The substrate polar- 
ization induced renormahzation of energy levels 
reduces the bandgap of substrate-supported gra- 
phone compared to the free-standing graphone. 

phone LDA bands in the vicinity of the valence- 
band maximum and the conduction-band mini- 
mum remain unaffected by the presence of hBN 
substrate. The shape of GW-corrected bands 
remains more or less unaffected, however, the 
bandgap of chair and boat-graphone when sup- 
ported on hBN substrate reduces by 0.86 eV and 
0.72 eV respectively compared to the bandgap of 
their free-standing counterparts. 

The reduction in GW bandgap is attributed to 
the polarization effects at the graphone/hBN in- 
terface. Similar reductions have been found in 
experiments and GW calculations for several dif- 
ferent molecules adsorbed on metals, semicon- 
ductors, and insulators. ^^"^^ Careful inspection of 
Figures [figure] [2] []2 and [figure] [3] []3 indicates 
that the GW corrections to the LDA bands are 
smaller in graphone on hBN compared to those in 
free standing graphone. This is because of the fact 
that the polarization of hBN substrate reduces the 
screened Coulomb potential, W, which in turn re- 
duces the GW bandgap. Similar bandgap reduc- 
tion due to the polarization of hBN substrate is 



expected in hydrogenated graphene with different 
hydrogen coverage. The bandgaps of GNRs are 
also reduced when they are supported on hBN sub- 
strate.^^ Non-local polarization effects can not be 
modeled in the LDA. The GW approach includes 
non-local polarization effects through the screened 
Coulomb interaction, however, it tends to be com- 
putationally extremely demanding. 

The bandgap reduction due to hBN substrate 
polarization can be estimated by using com- 
putationally much less demanding image-charge 
model. In this model, the quasiparticle en- 
ergy of a substrate-supported layer is given by 
E^^ , , = E^V + APj where E^f is the 

j;supported j;free ' J ^^^^'^^'^ j-Jree 

quasiparticle GW energy of state \j) of free stand- 
ing layer, graphone in this case, and APy is the cor- 
rection due to substrate polarization. The quasi- 
particle bandgaps of chair and boat-graphone esti- 
mated using the image-charge model are 2.03 eV 
and 4.36 eV respectively, which compare well 
with the bandgaps obtained using computationally 
extensive GW calculations including hBN sub- 
strate (Figure [figure] [3] [] 3). The details of image- 
charge model calculations and underlying assump- 
tions are discussed in the supplemental material. 

In addition to the bandgaps, the band offsets 
are of crucial importance for the successful use 
of a heterostructure in the electronic devices. As 
an illustration, we analyse the conduction and 
valence band offsets, denoted by AE'^ and IsEy, 
respectively, at the chair-graphone/hBN interface 
when hBN is placed on both the hydrogenated 
and non-hydrogenated sides of graphone. Fig- 
ure [figure] [4] []4(b) shows a schematic depict- 
ing the calculated band offsets in the hBN/chair- 
graphone/hBN heterostructure. Three monolayers 
of hBN are included on either sides of graphone. 
The band offsets do not change when more than 3 
monolayers of hBN are included. 

The bandgap of chair-graphone on hBN reduces 
from 1.93 eV (Figure [figure] [3] [] 3) to 1.74 eV in 
the presence of hBN on the hydrogenated side of 
graphone. This is due to the additional reduction 
of screened Coulomb interaction induced by the 
polarization of hBN layers on the hydrogenated 
side. The band offsets at chair-graphone/hBN 
interfaces are asymmetric such that the bands 
of hBN on the hydrogenated side of graphone 
are lower in energy compared to the bands of 
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Figure 4: Band offsets of chair- graphone em- 
bedded in hBN layers, (a) A supercell contain- 
ing 3 monolayers of hBN on both sides of chair- 
graphone. The direction of electric field induced 
by the dipoles in graphone layer is indicated by an 
arrow, (b) Band offsets at the graphone/hBN inter- 
faces. 



hBN on the non-hydrogenated side. To investi- 
gate the origin of this asymmetry, we carried out 
Bader charge analysis of a free standing chair- 
graphone. In chair-graphone, the hydrogenated 
C-atoms aquire a slight net positive charge while 
the non-hydrogenated C-atoms aquire a slight net 
negative charge creating a dipole layer with an 
electric field depicted by a thick grey arrow in 
Figure [figure] [4] [] 4(a). This electric field causes 
asymmetry in the band offsets at graphone/hBN 
interfaces. Similar modifications of the band off- 
sets have been found in conventional semiconduc- 
tor homo- and heteroj unctions when a dipole layer 
is inserted at the interface. 

In field effect devices employing chair-graphone 
as channel material, hBN can be used as a sub- 
strate as well as a dielectric layer separating the 
gate electrode and the graphone channel. The band 
diagrams in Figure [figure] [4] [] 4(b) indicate that 
hBN can be placed on either side of graphone 
when used as a substrate. When hBN is used as 
a dielectric layer, however, it should be placed 
on the hydrogenated side and not on the non- 
hydrogenated side of graphone. This is because 
of the fact that when hBN is placed on the non- 
hydrogenated side of graphone, AEy at the inter- 
face is too low (< 1 eV) to prevent the thermionic 
emission of holes into the hBN dielectric layer, 
which consequently results in gate leakage current 



degrading the performance of the field-effect de- 
vice. On the contrary, when hBN is placed on 
the hydrogenated side of graphone, it has suffi- 
ciently high barriers for both electrons (AEc) and 
holes (A£'v) to prevent their thermionic emission 
into the hBN dielectric layer. 

In summary, our first-principles calculations 
suggest that in a Bernal stacked graphene on 
hBN substrate a bandgap on the order of 0.1 eV 
opens up due to the sublattice asymmetry in 
graphene induced by hBN. This bandgap closes 
in the presence of slight misorientation from the 
Bernal stacking. Hydrogenation of graphene pro- 
vides a promising approach to open up bandgaps 
larger than 2.5 eV. The polarization effects due 
to surrounding hBN dielectrics, however, reduce 
bandgaps of graphone by about 1 eV. Thus, accu- 
rate electronic structure calculations of the atom- 
ically thin materials such as graphone should al- 
ways take into account the surrounding materi- 
als. The calculated band offsets suggest that in 
the field effect devices employing chair-graphone 
as a channel material, hBN can be used as a sub- 
strate as well as a dielectric layer separating the 
graphone channel and the gate electrode. Bandgap 
engineering of graphene by chemical functional- 
ization is currently an extremely active area of re- 
search and this work provides theoretical instruc- 
tion for analyzing the experimental observations 
such as the electronic structure modulation by the 
substrate. 
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